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Cosmic ray models compared to Fermi-LAT positron and electron separate 

spectra 
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In this symposium tiie Fermi-LAT collaboration released the results of its measurement of the e~ 
and separate spectra, obtained using the Earth magnetic field. Those results confirm PAMELA 
finding of an increasing positron fraction above 10 GeV and allow, for the first time independently 
on the electron spectrum, to reject cosmic ray spallation onto the ISM as the dominant positron 
production mechanism above that energy. We show, for two different propagation setups, that 
double component models which were developed to provide a consistent interpretation of PAMELA 
positron fraction and Fermi-LAT e"*" -|- e~ data, naturally reproduce also the Fermi-LAT e~ and 
separate spectra. 



I. INTRODUCTION 

One of the most striking recent results in cosmic 
rays (CR) physics was the observation performed by 
the PAMELA satellite experiment that the positron 
to electron fraction e+/(e^ + rises with energy 
from 10 up to 100 GeV at least [if. 

In the standard CR scenario, the bulk of electrons 
reaching the Earth in the GeV - TeV energy range 
are originated by Supernova Remnants (SNRs) and 
only a subdominant component of secondary positrons 
and electrons comes from the interaction of CR nuclei 
with the interstellar medium (ISM). Since the resi- 
dence time of primary nuclei in the Galaxy decreases 
with energy, in that scenario a growing positron frac- 
tion could be justified only by a very steep primary 
electron spectrum. This possibility, however, was 
ruled out by the Fermi-LAT collaboration which found 
the e~ -\- spectrum in the 7 GeV - 1 TeV energy 
range to be compatible with a power-law with index 
7(e±) = -3.08 ± 0.05 0, i], which is significantly 
harder than what estimated on the basis of previous 
measurements. Although excluded by the PAMELA 
collaboration, in principle, the rising positron frac- 
tion could be the consequence of misidentified CR pro- 
tons (which are 10"* 10^ times more abundant than 
positrons). Therefore, an independent confirmation of 
PAMELA results was called for. 

The awaited confirmation arrived just during this 
symposium as the Fermi-LAT collaboration released 
the preliminary results of the measurements of the 
absolute e+ and e~ spectra, and of their fraction, be- 
tween 20 and 120 GeV performed using the Earths 
magnetic field (W. Mitthumsiri et al. plenary talk). 
More recently, the collaboration published the final 
results of those measurements which were extended 
up to 200 GeV [3]. The sum of the e+ -I- e~ separate 
spectra was found to agree with the e+ -I- e~ spec- 
trum previously measured by Fermi-LAT. A steady 
rising of the positron fraction was observed by this 
experiment up to that energy in agreement with that 
found by PAMELA. In the same energy range, the 
e~ spectrum was fitted with a power-law with index 



7(e~) = — 3.19±0.07 which is in agreement with what 
recently measured by PAMELA between 1 and 625 
GeV Most importantly, Fermi-LAT measured, for 
the first time, the e"*" spectrum in the 20 - 200 GeV 
energy interval and showed it is fitted by a power-law 
with index 7(e+) = -2.77 ± 0.14. 

Remarkably, the e+ spectral index measured by 
Fermi-LAT is almost coincident with that of CR pro- 
tons which, in the standard scenario, are the main 
positron primaries. Since Feynman scaling implies 
that secondary positrons are produced with the same 
spectral index of protons, no room is then left to the 
steepening that propagation and energy looses should 
unavoidably produce in the propagated e"*" spectrum. 
From this simple argument it is already evident that 
the standard production scenario must be incom- 
plete. In the following we will explicitly show this is 
the case for two different propagation setups. We will 
then show to which extent and under which conditions 
the alternative scenario, which invokes the presence 
of a new primary positron and electron component, 
consistently describes all Fermi-LAT and PAMELA 
electron and positron data. 



II. SINGLE VS DOUBLE COMPONENT 
MODELS 

In 0, 0, 0| the Fermi-LAT collaboration discussed 
two scenarios to describe the e+-|-e~ spectrum it mea- 
sured. Schematically, we can define them as single and 
double component scenarios. 

The former scenario assumes one single class of 
sources, spatially distributed like Galactic SNR's, in- 
jecting e~ in the ISM with a, a priori unknown, bro- 
ken power-law spectrum. In those papers propagation 
was modeled using the GALPROP numerical diffusion 
package and solar modulation was then taken into ac- 
count in order to reproduce low energy data. This 
was done in the charge independent force field approx- 
imation by fixing the modulation potential $ against 
proton data taken in the same solar phase. 

Both models considered in [s*] assume a cylindri- 
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FIG. 1: The positron fraction computed with DRAGON or GALPROP is compared to experimental data. The shaded region 
represents the sum in quadrature of statistical and systematical Fermi- LAT errors. Upper/lower dotted lines represent the 
predictions of the standard single component scenario for strong/moderate reacceleration propagation setups described 
in Seclll fiinl Dashed/solid lines have been derived for double component models using the same propagation setups 
respectively. All models assume solar modulation with potential <& = 550 MV. 



cal diffusive halo with half-thikness of 4 kpc, a dif- 
fusion coefhcient scaling with rigidity like p^/^ (Kol- 
mogorov like diffusion) and relatively strong reacceler- 
ation (the Alfven velocity was taken va — 30 kms~^). 
Under those conditions, GALPROP provides an excel- 
lent description of most CR nuclei measurements @. 
It should be taken in mind, however, that the ob- 
served antiproton spectrum is not well fitted under 
those conditions while other choices of the propaga- 
tion parameters were shown to consistently reproduce 
both nuclear and antiprotons data (see Sec. IIIip . 

The single component reference model in Q is char- 
acterized by an e~ injection spectral index 70 (e~) = 
— 1.6/ — 2.5 below/above 4 GeV. That spectral break is 
required to reproduce low energy AMS-01 e~ data 10] 
as well as the spectrum of the synchrotron emission 
of the Galaxy below 1 Ghz [3. This model differs 
from similar pre-Fermi GALPROP models just for the 
harder spectrum it adopts above 4 GeV as required 
to track the Fermi-LAT + hard spectrum. An 
high energy cutoff in the source spectrum was also in- 
troduced in order not to overshoot H.E.S.S. e+ + 
data [l^l in the TeV region. In spite of such tuning 
this model does not allow a very satisfactory fit of 
the e"*" + e~ observed spectrum. Furthermore, since it 
assumes only e"*" secondary production, it cannot ex- 
plain the positron fraction rise observed by PAMELA 
and now confirmed by Fermi-LAT (see dotted lines in 
Fig. [T]). In Fig. [ifrsi we show that model also fails to 
describe the e+ measured by Fermi-LAT. This is the 
case also for the moderate reacceleration propagation 
setup that will be discussed in the next section. 

The double component scenario is characterized by 
the presence of a new term in the electron and positron 



source spectrum with the form 

Jcxtra(e±) cx E ^«(^*) exp(-£;/£;eut) . (1) 

For 7o(e^) ~ —1.5 and i?cut — 1 TeV this term has 
been shown to account not only for thepositron frac- 
tion anomaly observed by PAMELA Q but also an 
excellent fit of the -I- e~ spectrum measured by 
Fermi-LAT and H.E.S.S. [1, Q. Several hypothesis 
have been risen for its origin including acceler- 
ation in pulsar wind nebulae, dark matter annihila- 
tion, secondary production in SNRs (see [7j and 
Ref.s therein). Although the spatial distribution of 
the extra component source term generally depends 
on which of those scenarios is adopted, this has no 
consequences below few hundred GeV since at those 
energies the propagation length is comparable to 
the Galaxy size so that spatial features in the source 
term are aver aged-out. For this reason the following 
considerations apply both to astrophysical and to dark 
matter double component models. 

The double component model considered in Q 
adopts 7o(e~) = —1.6/ — 2.7 below/above 4 GeV for 
the standard electron component, 7o(e^) = —1.5 and 
-E-cut = 1.4 TeV for the electron and positron extra 
component. In Fig.s [T] and [3] we compare the pre- 
dictions of that model also with the new Fermi-LAT 
data ,4]. The reader can see from those figures as 
this model correctly reproduces not only the -t- e~ 
spectrum measured by Fermi-LAT but also its new e~ 
spectrum. The model also matches most Fermi-LAT 
e"*" data points and, most importantly, it reproduces 
their slope. The new Fermi-LAT data, therefore, con- 
firm the presence of an spectral component in the 
form given in Eq|T] 
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FIG. 2: The e"*" + e~ (black lines), the (red lines) and the e~ (blue lines) spectra computed for the single component 
scenario are compared with Fermi-LAT data in 3 • Dashed/solid lines have been derived using the strong/moderate 
reacceleration setups described in SeciII )flIII Statistical and systematic errors are summed in quadrature. All lines 
accounts for solar modulation with <I> = 550 MV. 
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FIG. 3: The predictions of the double component models discussed in Sees |lT] and IIIII are compared with Fermi-LAT 
data in 4] . The modulation potential and the line notation are the same of Fig. [21 



III. THE CASE OF MODERATE 
REACCELERATION 



In this section we test to which extent the choice of 
the propagation setup may affect our previous results. 
Our attempt is also motivated by the well known dis- 
crepancy between the prediction of strong reaccela- 
tion models and the positron fraction measured by 
PAMELA below 10 GeV. Such low energy PAMELA 
anomaly is commonly ascribed to solar modulation be- 
having in a more involved way than in the force field 
approximation. Charge dependent effects correlated 
to the solar magnetic field polarity are indeed theoreti- 
cally expected as a consequence of the complex spatial 
and temporal structure of the heliosphere (see [lal for 
a recent review). Typically, however, those effects are 
relevant only in the sub- GeV region and a satisfactory 



description of all available electron and positron data 
in those terms has not been yet achieved. 

Here we suggest that at least part of the low energy 
positron fraction anomaly could be related to other 
issues encountered by strong reacceleration propaga- 
tion setups: 1) as we mentioned, those models do not 
provide a satisfactory description of the antiproton 
spectrum measured by PAMELA and other experi- 
ments d, |T3|; 2) they reproduce the observed proton 
spectrum only by introducing a break ad hoc in the 
source spectrum; 3) strong reacceleration rise similar 
problems with the electron spectrum in the GeV re- 
gion [l^ ; 4) they are incompatible with the observed 
spectrum of the synchrotron emission of the Galaxy 
14| . Although none of these warnings might exclude 
strong reacceleration by itself, taken together they se- 
riously justify to consider alternative possibilities. 
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For all these reasons we consider here also a propa- 
gation setup characterized by moderate reacceleration 
{va = 15 kms"^) and a diffusion coefficient with a 
Kraichnan like dependence on rigidity: D{p) oc p^/^. 
More specifically we adopt here the KRA model pro- 
posed in p^ . That model was shown to correctly 
reproduce the B/C and other nuclear CR data if 
the particle speed dependence of the diffusion coeffi- 
cient is suitably tuned (as could be justified by MHD 
wave dissipation, see [9] and ref.s therein). Differently 
from strong reacceleration setups, it reproduces the 
antiproton spectra measured by PAMELA and, 
what is most relevant here, it requires no break in the 
proton source spectrum at least up to few hundred 
GeV. In fact, we verified that this is crucial to cor- 
rectly match the PAMELA positron fraction below 
10 GeV. Under those conditions, and assuming the 
presence of a extra-component as that adopted 
in the previous section, the + e~ spectrum mea- 
sured by Fermi-LAT, H.E.S.S. and AMS-01 was also 
accurately reproduced. Respect to Ref. [lB| here we 
only change the low energy e~ source spectral index 
by taking 7o(e^) = —1.6/ — 2.65 below/above 4 GeV 
rather that —2/2.65. In Fig. [l]the reader can see as 
this arrangement allows not only a good fit of the e"*", 
e~ and e+/(e+-|-e~) Fermi-LAT data between 20 and 
200 GeV but also to improve considerably the descrip- 
tion of the PAMELA positron fraction data below 20 
GeV. The small residual offset may either be due to 
solar modulation or to a experimental systematic (the 
reader may have noted the small offset between the 
6+ fraction measured by PAMELA and Fermi-LAT). 

We verified that PAMELA e~ data [5j are also con- 
sistently reproduced by the model considered in this 
section, with almost the same modulation potential, 
if their normalization is rescaled so to make them to 
coincide with Fermi-LAT e~ data at 20 GeV. 



IV. CONCLUSIONS 

In this symposium the Fermi-LAT collaboration 
presented its measurements of the positrons fraction 



as well as of the e+ and absolute spectra in the 
20 - 200 GeV energy interval. We showed that the 
spectrum measured by Fermi-LAT is in contrast 
with the standard scenario in which are originated 
only by the spallation of the nuclear component of CR 
onto the ISM. This evidence is even stronger than that 
provided by the positron fraction anomaly found by 
PAMELA, and confirmed by Fermi-LAT, since it is in- 
dependent on the e~ spectrum reducing the involved 
systematics. We then tested the alternative scenario 
which invoke the presence of an additional electron 
and e+ primary spectral component, against the same 
data. We showed, for two different propagation se- 
tups, that models which were proposed in this frame- 
work to consistently interpret PAMELA positron and 
the Fermi-LAT e"*" -I- e~ spectrum naturally reproduce 
the e"*" and e~ separate spectra. By measuring the 
e+ and e" spectra above 200 GeV, AMS-02 will be 
necessary to have a clue on the nature of the extra 
component discovered by PAMELA and Fermi-LAT. 
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